*  ARKIV FOR FYSIK Band 16 nr 30 
f Communicated 26 November 1958 by Ivan WALLER 


Cluster development of the density matrices in the Jastrow 


method for many-fermion systems 


By S. O. LunpDevist 


ABSTRACT 


A method is given for the cluster development of the reduced density matrices for Jastrow 
wave functions. The essence of the method is to make extensive use of the simple properties of 
the Dirac density matrices corresponding to the uncorrelated motion, in conjunction with the 
method of expanding the distribution functions for a classical imperfect gas given by Mayer and 
Montroll. 


Introduction 


In the original paper by Jastrow [1], he discusses the case of a uniform Fermi gas. 
Starting’ from the cluster expansion of the pair distribution function for a classical 
gas and taking the antisymmetry requirement into account in an approximate way, 
he derived an approximate expression for the energy of the gas. A more general 
derivation, not restricted to uniform systems and taking account of the exclusion 
principle, was later given by Iwamoto and Yamada [2], who made use of a trans- 
formation introduced by Ursell [3] in the theory of imperfect gases. They only con- 
sider the cluster development of the energy so no information about the distribution 
functions is obtained from their work. Recently a detailed discussion has been given 
by Hartogh and Tolhoek [4] of the distribution functions for a uniform Fermi gas, 
neglecting spin. 

In the present note a different method will be indicated, which is not restricted 
to a uniform gas and in which spin and isobaric spin is taken into account from the 
beginning. The method has some similarity to that used by Hartogh and Tolhoek, 
but differs from theirs by making extensive use of the simple properties of the Dirac 
density matrices for the system [5]. No special attention has to be given to the fulfil- 
ment of, e.g., conservation of linear or angular momentum components, because such 
conservation principles are already contained in the Dirac matrices, if the one- 
particle functions are properly chosen. 


The reduced density matrices for Jastrow wave functions 


We consider a system of A fermions with equal masses m, which interact with 
two-body forces described by a potential v = v(i,j), where i and j stand for the 
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coordinates (including spin and isospin) of the particles 7 and j. Let us introduce 
the reduced density matrix for n particles [6] 


(1’2’...n'|o|12...0)= 
A! PES 20. 50, ward a AYP (i @E1 2 Alar 


= (1) 


(A—n)! Pe (ea lee 


where y(1....A) is the wave function for the system. dr implies summation over 
discrete as well as integration over continuous variables. The normalization of (1) 
is chosen such that the diagonal element represents the probability density for n 
particles. The expectation value of the energy can be expressed by means of the 
reduced density matrices for one and two particles as follows: 


f 1 
B= 5 |dt,(1|pie|1)+5 Jdt,dr_(12|v9|12). (2) 


Here p, denotes the momentum operator of particle one and the symbols (1 | pie | 1) 
and. (12|ve|12) mean ordinary matrix products. 
As a trial function for the ground state of the system we write 


p(1... A) =F O(1... A), (3) 


where ®(1... A) describes the uncorrelated motion and is taken as a Slater deter- 
minant of suitably chosen one-particle functions ¢, (2), and F is an operator describing 
the correlations. 

In the Jastrow method one chooses 


7 pails. (79), (4) 
where the two-particle correlation function f(ij) satisfies the symmetry condition 
f(%j) =f (92) (5) 

and the boundary condition 
f(j)—1, ty —-Co. (6) 


We shall restrict ourselves to the case where [(ij) depends only on the position and 
spin and isospin of the particles. It would be desirable to include also a dependence 
on the states of the particles 7 and j. In such a case (4) has to be replaced by the 
symmetrised product of the f(ij). This more general case can easily be treated by 
an extension of the method used here. 


In the following we shall make use also of the density matrices, corresponding to the 
uncorrelated motion [5] 


(Urn lo|...2>= 


A! 
Wwe for 6 WED ADOUL PnP) od)dr4 (7) 


if ® is normalized to unity. 
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xpansic of (1... n sth 1...) can be performed by a simple extension of the 
S erbdubed by Mayer and Montroll [7] for an eee classical gas. It is 

is snvenient to introduce the notation 
: —Q(A-1= TT 99), (11) 
_ where Sandie g (*9) =/* (04) f(t) =h(aj) + 1. (11’) 


_ To simplify the notation the primes on certain quantities will not be written out in 


the following. For instance, if i <j <1, g(t) means f* (i’9’) (17) and if 1 <1 <j, g (7) 


means f*(i’j)f (27). Under the assumption of short range Ned the quantities h (77) 


ae a 


tend rapidly to zero with increasing 74; 


Writing F* F=Q(A)= zal g (i) I (1+A(kl))Q(A—n) (12) 
os ees 


and expanding I1(1 + h(kl))Q(A —n) we get an expansion of the form 


_s_(4=")! am es 
QA)=S Gaza HP wD O(n —A), (13) 
- where 
GO (1... n+ A)= GS (1... m) ETA (ES) 13’) 
Syith OP (L..m)= TT gli) n>1 
and GP (1... n=l. 


* 


" Introducing the expansion (13) into (10) we get 


-. 


| 


ate. . Ee a , 2 Shy P 

(A—n)! | | , aie 

wis DIAS i) he a Al epee ra Be 
me JQ(A)... Afg|1... Aydr4 


(U2. Wott 25g m= 2. 


(14) 


Because of the effect of exchange between the particles 1...» +A and the remain- 
ing particles, we cannot factorize (14) directly as in classical theory. Instead we can 
proceed as follows by successive expansions. For the sake of brevity we first introduce 
the shorter notations 


On =(1’...n'lo|1...2) and o7 =<1'...n'lo|1...m>. 


I. In the first step we make use of the Laplace expansion of a determinant, which 
applied to 04 =det,<x;|o|2;> gives 


04 = 0n+104-n-a+ Raa, . (15) 


where R,, is the sum of all the terms representing exchanges between the particles 
1...n +A and the remaining ones. The first term permits a factorisation and gives the 
contribution 


DZ(n+)[ QP (L...ntAorude (16) 
4a 


A= ee | A-—m 
with Y, (m) = _ At [Q(A=moh-mdrt™ . (16’) 


(A—m)! ff Q(A) oR dr4 


We note that if we expand Q(A) in the denominator in (16’) and apply (15) with 
neglect of R,, we get for n =2 a trivial generalisation of Jastrow’s derivation of 
the pair distribution function. In this approximation all exchanges between the 
particles 1.... +A are taken into account, but the exchanges with the remaining 
particles are neglected. These exchanges give exchange couplings between clusters, 
which are disconnected in classical theory. Such effects are taken into account by 
Iwamoto and Yamada [2] and have recently been discussed by Hartogh and Tolhoek 
[4]. The exchange effects are contained in the terms R,, in (15). The Laplace expan- 
sion is of no direct use as Q(A —n— A) contains all kinds of clustering between the 
particles, so we use instead the following procedure. 
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over the factor Q(A —n—A-—w). We have then, however, included also terms which 


7?’ 


7 rs 
fie 


ae | 
the index 0 indicating that there is no connexion with the basis 1.... We wish to 


peal) hh 13) (hy Te ee 


(A=n—A)! 


Soa CY Res a 


GO (ntA+1...mtAtp)Q(A—n—-A—p) 


consider exchanges between the particles 1...n +A andn+A+1...1+A+4 only. 
Then applying the Laplace expansion as in I we can for each A separate the integral 


only give interchanges within the two sets of coordinates. Those terms are easily 


oy 


eliminated by replacing gf 1+. by 


(olan = Ontie y= CFs 108: (18) 


For each A we then get the contribution 


=1 


FS (n+ A+ py) [ OP (L... mt AGP (mt At]... n+A+p)[Olnsandt™. (19) 
i : 


III. These arguments can now be repeated to include exchange coupling between 
an arbitrary number of disconnected sets of particles. To consider exchange coupling 
between three disconnected sets we replace for each 4 and wu Q(A —n—A-—w) in 
(17) by 


pF AH EAL MEAT BEN QA mA Hm9) (20) 
y=1 Rie Cele Ee FILL: : > 


Factorising the expression as in II we get for each A and y the three-clusters exchange 
contributions expressed as 


DZ(ntA+ ptr) [ GPL... mt AGP (ntAtl..ntatp)x 
v=1 


XG (nt At mt l...mtAt+et+)[Olnsrirnrvd tl”, (21) 


Here [0]n+2++» is determined so that in each diagram the three sets become connected 
if we represent each exchange by a certain type of link [4]. The expression for [0]n+2++» 
is easily written down by subtracting from 0n+1++» those components which corres- 
pond to forbidden diagrams. Denoting the three sets by 1, 2 and 3 for simplicity, 
we get the following forbidden diagrams. 
1) 1,2 and 8 all disconnected, which corresponds to gf @2’@3 (no exchange between 
the groups). 
2) Two sets (e.g. 1 and 2) connected by exchange, and the third (3) disconnected, 
which corresponds to 03[@1+2 —@102] and correspondingly for the other two 
cases. 
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ise be 


A=0 ps 


4 ; : ; ae a Sean age (0) 7 use : 
+S DD Avtar uty) [GP mt AEN (ntAt1 . 
BO ast ‘ » +2 ; 4i® e rah: % 
XGP (nt dt wtl...mth+ utr) [Olnvarured MY + veaioad wae 


; i 
This formula can now be simplified essentially by taking together certain classes 
of terms. We note that in deriving (22) we make use of formulas like (17) and (19) — 
where the terms w =1, »=1,... represent a cluster consisting of one single point. 

Such terms in (22) which contain single points can be reduced noting that the inte- 
grands contain factors of the type stat. 


<x;|0 | tq><te|0 |x, 


x, being the coordinate for an isolated point. As the integrand does not depend on 
x, in any other way, the integration over x, gives the result 


<xi|e|x)> 


which means that the integral reduces to one where this particular one-particle 
cluster does not appear. In this way every integral with a number of single-point 
clusters can be reduced (except for a numerical factor) to a term in (22) with no 
single-point clusters. Thus every term in (18) without any single-point clusters can 
be regarded as the parent of a family of terms obtained by attaching exchange links 
to a number of single points. The contribution from each member of the family 
differs only by a factor from that of the parent so the sum over the whole family 
will be factor times the contribution from the parent. Summing up over all such 
families we still have an expansion of the form (22) but with the factors Z(m) 
replaced by new factors Z’(m) and without the terms u =1, »=1... 


Qn = 2’ (n) G9 (L...m)or + > 2! (n+) [ QPL... w+ A) oPeadt?+ 


+ 2, 2, A (m+ At py) f OP (1... m+ AGP (m+ 441... 41+ w) [ola rrend e+ 


+> > ; Z' (n+ At pwtv) { QP (L...mt4GO(ntat1...ntAtp)x 


XG (M+ At WAL... W+A+ U+9)[Olnsatury DH? +. (23) 
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s terms which tend to zero when to infinity, ¢ 
Mee et Av 
2'(m= [1+ 5 aimed ey 


The constants «{” depend on the strength of the two-body correlation through a 

‘parameter, which for a uniform gas is simply 7 = Aw,/Q, where w, is the weighted 
_ volume in which the correlation is appreciable and Q is the volume. «{” tend to 
_ zero when 7 goes to zero. Further for m< A it is easily found that 


4 Z' (m+1) ‘ 
ce a a Ge , 
— Z’ (m) (independent of m) (24’) 


_ Assuming that the expansions converge rapidly enough only terms fulfilling this 
condition will be of importance and then only powers of Z’ will appear in (23). 
The formulas given above are obtained by first taking the single-point clusters into 
_ account and then reducing the formulas. If we simply had left out these single-point 
_ clusters from the beginning, we would have obtained (23) but with the original 
_ normalisation constants Z(m) instead of Z’(m). It is now simple to show that if we 
treat the denominator in (16’) in the same way, i.e. neglect all single points in the 
expansion of Q(A), the original Z(m) satisfies the relation (24) with the same «{” 
~ and further that Z(m + 1)/Z(m) =Z (independent of m). From this we can conclude 
that Z’ = Z, i.e. (23) is just the expansion we would have obtained if all exchanges 
involving point clusters had just been left out from the very beginning. This remark 
- connects our approach with that by Hartogh and Tolhoek for a uniform gas, and 
especially we can apply their result that Z—1 as A—oo. 
The final expression will thus take the form 


On= 2G (1 Mon + >, ZP*? | GSP (1... m+ A)onsadt' + 
=1 
+> > Zrte GPL... n4AGD (Wt Atl... nt+At p)[olnzarnde ? + 
A=0 p=2 


ee Aa GOL 3. 0A) GD (ot AL... w+ Ap) X 


A=0 M=2 v=2 
x GO (ntAt+prl... rhe er?) Cie dto ven (25) 


It is now a trivial matter to write down the first few terms in the expansions of 
0, and g2, from which one then gets the expansion of the energy. We do not give 
those rather lengthy formulas here, as we have here wished to indicate this method 
of deriving the expansion only, rather than give ali the details. More of the details 
will be given in a forthcoming paper. 

Order of magnitude estimates of the terms can easily be made for a uniform 
system by noting that the one-particle Dirac matrix <a'|o|a> isa smooth decreasing 
function of |r’—r| proportional to A/ Q and that all exchange terms give contribu- 
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